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ABSTRACT 

The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neigh- 
borhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic 
capabilities of the mission are being used to determine precise radii and ages for the target stars from their 
solar-like oscillations. Chaplin et al. (2010) published observations of three bright G-type stars, which were 
monitored during the first 33.5 d of science operations. One of these stars, the subgiant KIC 11026764, ex- 
hibits a characteristic pattern of oscillation frequencies suggesting that it has evolved significantly. We have 
derived asteroseismic estimates of the properties of KIC 11026764 from Kepler photometry combined with 
ground-based spectroscopic data. We present the results of detailed modeling for this star, employing a vari- 
ety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints 
simultaneously. We determine both the radius and the age of KIC 1 1026764 with a precision near 1%, and an 
accuracy near 2% for the radius and 15% for the age. Continued observations of this star promise to reveal 
additional oscillation frequencies that will further improve the determination of its fundamental properties. 

Subject headings: stars: evolution — stars: individual (KIC 11026764) — stars: interiors — stars: oscillations 
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1. INTRODUCTION 

In March 2009 NASA launched the Kepler satellite, a mis- 
sion designed to discover habitable Earth-like planets around 
distant Sun-like stars. The satellite consists of a 0.95-m tele- 
scope with an array of digital cameras that will monitor the 
brightness of more than 150,000 solar-type sta rs with a pre- 
cision of a few parts-per-million for 4-6 years (Borucki et al. 
l2007h . Some of these stars are expected to have planetary sys- 
tems, and some of the planets will have orbits such that they 
periodically pass in front of the host star, causing a brief de- 
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Figure 1. Evolution of the / = (dotted) and 1 = 1 (solid) oscillation frequencies as a function of age for a representative stellar model of KIC 1 1026764. The 
frequency separation between consecutive / = 1 modes during an avoided crossing is a strong function of the stellar age. Note the prediction of a high-frequency 
avoided crossing above 1250 fiYlz for Model AA, indicated by the vertical line. 



crease in the amount of light recorded by the satellite. The 
depth of such a transit contains information about the size of 
the planet relative to the size of the host star. 

Since we do not generally know the precise size of the 
host star, the mission design includes a revolving selec- 
tion of 512 stars monitored with the higher cadence that 
is necessary to detect short period solar-like oscillations, 
allowing us to apply the techniques of asteroseis mology 
dChristensen-Dalsgaard etaT]|2007UAerts et alj|2010l) . Even 
a relatively crude analysis of such measurements can lead to 
reliable determinations of stellar radii to help characterize the 
extra-solar planetary systems discovered by the mission, and 
stellar ages to reveal how such systems evolve over time. For 
the asteroseismic targets that do not contain planetary com- 
panions, these data will allow a uniform determination of the 
physical properties of hundreds of solar-type stars, providing 
a new window on stellar structure and evolution. 

Initial results fr om the Kepler Asterose ismic Investigation 
were presented by Gillil and et al.l (l2010al) . while a more de- 
tailed analysis of the solar-like oscillations detected in several 
early targets was published by Ch aplin et alj d2010t) . The lat- 
ter paper includes observations of three bright (V~9) G-type 
stars, which were monitored during the first 33.5 d of science 
operations. One of these stars, the subgiant KIC 1 1026764 (= 
2MASS J19212465+4830532 = BD+48 2882 = TYC 3547- 
12-1), exhibits a characteristic pattern of oscillation frequen- 
cies suggesting that it has evolved significantly. 

In unevolved stars, the high radial order (n) acoustic os- 
cillation modes (p-modes) with a given spherical degree (/) 
are almost evenly spaced in frequency. As the star evolves 
and the envelope expands, the p-mode frequencies gradu- 
ally decrease. Meanwhile, as the star becomes more cen- 
trally condensed, the buoyancy-driven (g-mode) oscillations 
in the core shift to higher frequencies. This eventually leads 



to a range of frequencies where the nonradial (I > 0) oscil- 
lation modes can take on a mixed character, behaving like 
g-modes in the core and p-modes in the envelope ("mixed 
modes"), with their freq uencies shifted as they undergo so - 
called avoided crossings (lOsakilll975HAizenman et1d1ll977h . 
This behavior changes relatively quickly with stellar age, 
and propagates from one radial order to the next as the star 
continues to evolve (see Figure [T). Consequently, those 
modes that deviate significantly from uniform frequency spac- 
ing yield a strong (though m odel-dependent) constraint o n 
the age of the star (e.g., see Christensen-Dalsgaard 2004). 
Avoide d crossings have been observed in the subgian t stars 
j] Boo dKieldsenetal.1119951. 120031; ICarrieret alj 12005b and 
B Hvi (|Bedding et al.l 120071) . and p ossibly also in Procyo n 
(Bedd ing et alj|2010t) and HP 4938 5 (iDeheuvels et alJl2010t) . 
As noted bv iGilliland etall (12010a!) and lChaplin et al.l (12010) 
the dipole (/ = 1) modes observed in KIC 1 1026764 show the 
signature of an avoided crossing, raising the exciting possi- 
bility that detailed modeling of this star will provide a very 
precise determination of its age. 

In this paper we derive the stellar age, radius and other char- 
acteristics of KIC 11026764 by matching both the observed 
oscillation frequencies from Kepler photometry and the best 
available spectroscopic constraints from ground-based obser- 
vations. We describe the extraction and identification of the 
oscillation frequencies in §|2] and the analysis of ground-based 
data for spectroscopic constraints in $3] In $4] we provide the 
details of the independent codes and analysis methods used 
for the fitting, and in f5] we describe our final modeling re- 
sults. We summarize and discuss the broader significance of 
the results in §|6] 

2. OSCILLATION FREQUENCIES 

The 58.85-second (short-cadence) photometric data on 
KIC 11026764 came from the first 33.5 d of science opera- 
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Table 1 

The minimal and maximal sets of observed oscillation frequencies for KIC 1 1026764. 



* 

n 


Minimal Frequency Set (/iHz) 


Maximal Frequency Set (/iHz) 


/ = 


/ = 1 


/ = 1 


/ = 


7 1 
/ = 1 




10 












615.49 ±0.45 b 


11 


620.79 ±0.29 c 


653.80 ±0.23 




620.42 ±0.37 c 


654.16 ±0.39 


670.23 ±1.21 a 


12 








673.97 ±0.71 a 


699.87 ±0.66 b 


716.36±0.21 b 


13 


723.63±0.10 


755.28 ±0.26 


767 .46 ±0.46 


723.30 ±0.24 


754.85 ±0.23 


769.16±0.71 


14 


772.82 ±0.30 


799.96 ±0.40 


817.91 ±0.58 


772.53 ±0.33 


799.72 ±0.23 


818.81 ±0.30 


15 


822.72 ±0.14 


847.57 ±0.30 


867.66 ±0.86 


822.46 ±0.31 


846.88 ±0.34 


868.31 ±0.24 


16 


873.55±0.14 


893.48 ±0.33 




873.30 ±0.27 


893.52 ±0.20 


919.31 ±0.52 a 


17 


924.53 ±0.37 


953.57±0.39 


969.77 ±0.36 


924. 10 ±0.29 


953.51 ±0.22 


970.12 ±0.68 


18 


974.59 ±0.35 


1000.41 ±0.52 


1020.72 ±1.33 


974.36 ±0.26 


1000.38±0.41 


1019.72 ±0.44 


19 


1025.48 ±0.63 


1049.99 ±0.35 




1025.16±0.37 


1049.34 ±0.29 


1072.49 ±0.62 a 


20 


1076.70 ±0.29 






1076.52±0.51 







Reference value of n, not used for model-fitting. 
Mode not present in any of the optimal models. 



* Observed mode adopted for refined model-fitting. 
Models suggest an alternate mode identification (see ^6) 



tions (2009 May 12 to June 14). Time series data were then 
pr epared from the raw ob servations in the manner described 
by iGilliland et al.1 (1201 Obi). The powe r spectrum is shown in 
Figures 1 and 2 of Chapli n et al.l (120101) . Eight teams extracted 
estimates of the mode frequencies of the star. The teams used 
slightly different strategies to extract those es timates, but the 
main idea was to maximize the likelihood (Anders on et al.l 
1990) of a multi-parameter model designed to describe the 
frequency-power spectrum of the time series. The model in- 
cluded Lorentzian peaks to describe the p-modes, w ith flat and 
power-law terms in frequency (e.g., Harvey 1985) to describe 
instrumental and stellar background noise. 

The fitting strategies followed well-established recipes. 
Some teams performed a global fit — optimizing simultane- 
ously every free par ameter needed to descri be the observed 
spectrum (e.g., see Appourch aux et al.1 120081) — while others 
fit the spectrum a few modes at a time, an a pproach tradi- 
tional ly adopted for Sun-as-a-star data (e.g., see lChaplin et al.l 
1999). Some teams also incorporated a Bayesian approach, 
with the inclusion of priors in the optimization and Markov 
Chain Monte Carlo (MCMC) analysis to map the poste- 
rior distributions of the estimated frequ encies (e.g., see 
iBenomar et al.ll2009b ICampante et al.ll2010l) . 

We then implemented a procedure to select two of the eight 
sets of frequencies, which would subsequently be passed to 
the modeling teams. Use of individual sets — as opposed to 
some average frequency set — meant that the modeling could 
rely on an easily reproducible set of input frequencies, which 
would not be the case for an average set. We selected a mini- 
mal frequency set to represent the modes that all teams agreed 
upon within the errors, and a maximal frequency set, which 
included all possible frequencies identified by at least two of 
the teams, as explained below. 

From the sets of frequencies /v„/ , provided by the eight 
teams, we calculated a list of average frequencies v n \. For 
each mode {«,/}, we computed the number of teams return- 
ing frequencies that satisfied 



\v n lj-Vnl \ < PnlA 



(1) 



with <r„i j representing the frequency uncertainties returned by 
each team. We then compiled a minimal list of modes. For 



each {n,l} we counted the total number of teams with iden- 
tified frequencies, as well as the number of those frequen- 
cies that satisfied Eq.([T|i. Modes for which all identifications 
satisfied the inequality were added to the minimal list. We 
also compiled a maximal list of modes, subject to the much 
more relaxed criterion that the {«,/} satisfying Eq.([TJ should 
be identified by at least two teams. 

In the final stage of the procedure, we computed for each 
of the eight frequency sets the normalized root-mean-square 
(nns) deviations with respect to the v n \ of the minimal and 
maximal lists of modes. The frequency set with the smallest 
rms deviation with respect to the minimal list was chosen to 
be the minimal frequency set, while the set with the smallest 
rms deviation with respect to the maximal list was chosen to 
be the maxim al frequency se t . The minimal frequency set was 
also used by IChaplin et a l. (2010), and provided the initial 
constraints for the modeling teams (see Table [TJ. The max- 
imal frequency set was used later for additional validation, 
as explained in §0 Note that the same modes have slightly 
different frequencies in these two sets, since they come from 
individual analyses. The true radial order (n) of the modes 
can only be determined from a stellar model, so we provide 
arbitrary reference values for convenience. 

3. GROUND-BASED DATA 

KIC 11026764 (a 20 oo = 19 h 21 m 24?65, 5 2 ooo = 
+48°30'53."2) has a magnitude of V = 9.55. The atmo- 
spheric parameters given in the Kepler Input Catalog 36 
(KIC; lLatham et al.l 12005b as derived from photometric 
observations acquired in the Sloan filters are T & s = 5502 K, 
logg = 3.896 dex, and [Fe/H]= -0.255 dex. The quoted 
uncertainties on these values are 200 K in T e ff and 0.5 dex in 
\ogg and [Fe/H]. Since this level of precision is minimally 
useful for asteroseismic modeling, we acquired a high- 
resolution spectrum of the star to derive more accurate values 
of its effective temperature, surface gravity, and metallicity. 

3.1. Observations and Data Reduction 

The spectrum was acquired with the Fibre-fed Echelle 
Spectrograph (FIES) at the 2.56-m Nordic Optical Telescope 

36 http://archive.stsci.edu/kepler/kepler_fov/seai'ch.php 
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(NOT) on 2009 November 9 (HJD 2455145.3428). The 
1800 s exposure covers the wavelength range 3730-7360 A 
at a resolution R^67000 and signal-to-noise ratio S/N= 80 
at 4400 A. The Th-Ar reference spectrum was acquired im- 
mediately after the stellar spectrum. The reduction was per- 
formed with the FIESTOOL software, which was developed 
specifically for the FIES instrument and performs all of the 
conventional steps of echelle data reduction 37 . This includes 
the subtraction of bias frames, modeling and subtraction of 
scattered light, flat-field correction, extraction of the orders, 
normalization of the spectra (including fringe correction), and 
wavelength calibration. 

3.2. Atmospheric Parameters 

We derived the atmospheric parameters of KIC 11026764 
using several methods to provide an estimate of the exter- 
nal errors on 7^, logg, and [Fe/H], which would be used in 
the asteroseismic modeling. The five inde pendent reductions 
includ ed: the VWA 38 softwa re package dBruntt et all 12004 
I20l0ah. the MOOG 39 code dSnedenll 19731) . th e ARES 40 code 
(ISousaetalJl2007l). th e SYNSPEC method (lHubenvlll988l: 
Hubenv&Lanzlll995l) . and the ROTFIT code dFrasca et all 
2003112006b . The principal characteristics of the methods em- 
ployed by each of these codes are described below, and the 
individual results are listed in Table [2] 

For the VWA method, the T e s, \ogg and mi croturbulence of 
the ad opted MARCS atmospheric models (Gustafsso n et al.1 
2008) are adjusted to minimize the correlations of Fel with 
line strength and excitation potential. The atmospheric pa- 
rameters are then adjusted to ensure agreement between the 
mean abundances of Fe I and Fe II. Additional constraints on 
the surface gravity come from the two wide Ca lines at 6122 
and 6162 A, and from the Mg-lb lines dBruntt et al.ll2010bl) . 
The final value of \ogg is the weighted mean of the results 
obtained from these methods. The mean metallicity is cal- 
culated only from those elements (Si, Ti, Fe and Ni) exhibit- 
ing at least 10 lines in the observed spectrum. The uncertain- 
ties in the derived atmospheric parameters are de termined by 
perturb ing the computed models, as described in Bru ntt et al.l 
(2008). Having computed the mean atmospheric parameters 
for the star, VWA finally determines abundances for all of 
the elements contained in the spectrum (see Figure [2}. No 
trace of Li I 6707.8 absorption is seen in the spectrum of 
KIC 11026764. We estimate an upper limit for the equiva- 
lent width EW < 5 mA. 

The 2002 version of the MOOG code determines the iron 
abundance under the assumption of local thermodynamic 
equilibr i um (L TE), using a grid of ID model atmospheres by 
Kurucz (1993). The LTE iron abundance was derived from 
the equivalent widths of 65 Fe I and 10 Fe II lines in the 4830- 
6810 A range, measured with a Gaussian fitting procedure 
in the IRAF 41 task splot. For the an alysis, we followed the 
prescrip tion of | Randich et ail d2006l) . using the same list of 
lines as iBiazzo et al.l d2010l) . The effective temperature and 
microturbulent velocity were determined by requiring that the 

37 http://www.not.iac.es/instruments/fies/fiestool/FIEStool.html 

38 http://www.hans.bruntt.dk/vwa/ 

39 http://verdi.as.utexas.edu/ 

40 http://www.astro.up.pt/~sousasag/ares/ 

41 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



Table 2 

Atmospheric parameter estimates for KIC 1 1026764. 



r ef f logg [Fe/H] Method 

5640±80 3.84±0.10 +0.02±0.06 vwa 

5750±50 4.10±0.10 +0.11 ±0.06 MOOG 

5774 ±39 4.01 ±0.07 +0.09 ±0.03 ares 3 

5793 ±26 4.06 ±0.04 +0.10 ±0.02 ares' 3 

5630±70 3.79±0.17 +0.10±0.07 synspec 

5777 ±77 4. 19 ±0.16 +0.07 ±0.08 rotfit 



a 40 Fe I and 1 2 Fe II lines from lSousa eTal] <2T)06|) . 
b 247 Fe I and 34 Fe II lines from lSousa eTal] [2008). 



iron abundance be independent of the excitation potentials 
and the equivalent widths of Fel lines. The surface gravity 
was determined by requiring ionization equilibrium between 
Fe I and Fe II. The initial values for the effective temperature, 
surface gravity, and microturbulence were chosen to be solar 
(T eS = 5770 K, logg = 4.44 dex, and £ = 1 . 1 km s" 1 ). 

ARES provides an automated measurement of the equiv- 
alent widths of absorption lines in stellar spectra: the LTE 
abundance is determined differentially relative to the Sun with 
the help of MOOG and a grid o f ATLAS -9 plane-parallel 
model atmospheres jKuruczl|19 93). We used ARES with two 
different li sts of iron lines: a 'short' one composed of isolated 
iron lines dSousa et al.ll2006l) . and a 'long' one c omposed of 
iron l ines suitable for automatic measurements (Sou sa et all 
120081) . Our computations resulted in two consistent sets of 
atmospheric parameters for KIC 1 1026764. 

SYNSPEC provides synthetic spectra based on model at- 
mospheres, either calculated by TLUSTY or taken f rom the 
literature. We used the new gri d of ATLAS -9 models dKuruczl 
119931: ICastelli & Kuruczll2003l) to calculate synthetic spectra, 
which were then compared t o the observed spec trum. Based 
on the list of iron lines from lSousa et al.l d2008l) . we derived 
the stellar parameters in two ways to estimate the uncertainty 
due to the normalization of the observed spectrum. For the 
first approach we determined the minimum \~ of the deviation 
between the synthetic iron lines and the observed spectrum for 
a fixed set of stellar parameters. For the second approach, we 
determined the best fitting effective temperature and surface 
gravity for each iron line from the list, and adopted stellar pa- 
rameters from the mean. The primary uncertainty in the final 
parameters arises from the correlation between the effective 
temperature and surface gravity: a reduction of the effective 
temperature can be compensated by a reduction of the gravity. 
While the \ 2 method results in lower values for both parame- 
ters (T e ff = 5560 K, logg = 3.62 dex), the averaging approach 
yields higher values (r e ff = 5701 K, logg = 3.95 dex). The dif- 
ferences between the two results exceed the formal errors of 
each method. We therefore adopt the mean, and assume half 
of the difference for the uncertainty. The metallicity is de- 
termined by minimizing the scatter in the parameters derived 
from individual lines. 

ROTFIT performs a simultaneous and fast determination of 
r e ff, logg, and [Fe/H] for a star — as well as its projected ro- 
tational velocity vsin/ — by comparing the observe d spectrum 
with a library of spectra fo r reference stars (see iKatz et all 
119981: ISoubiran etail 119981) . The adopted estimates for the 
stellar parameters come from a weighted mean of the pa- 
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Figure 2. The mean abundances of the elements in the spectrum of KIC 1 1026764 as derived with the VWA software, including neutral lines (circles) and singly 
ionized lines (squares). The species labeled in bold are alpha eleme nts. The horizonta l bar indicates the mean metallicity of the star and the Icr uncertainty range 
of the determination. The abundances are given relative to solar ( Grevesse at al. 2007). 



rameters for the 10 reference stars that most closely resem- 
ble the target spectrum, which is quantified by a % 2 measure. 
We applied the ROTFIT code to all echelle orders between 
21-69, which cover the range 4320-6770 A in the observed 
spectrum. We also derived a projected rotational velocity for 
KIC 11026764 of 2.8 ±1.6 km s" 1 . 

3.3. Adopted Spectroscopic Constraints 

The effective temperature of KIC 11026764 derived from 
the five methods outlined above generally have overlapping 
lcr errors. They all point to a star that is hotter than the 
KIC estimate by 1 00-300 K. We find reasonable agreement 
between the derived values for surface gravity and the KIC 
estimate. Most of the applied methods result in \ogg above 
4.0 dex, again slightly higher than in the KIC. Only VWA 
and SYNSPEC yield slightly lower values, which is not sur- 
prising considering the correlation between T e ff and \ogg. Fi- 
nally, all of the methods agree that the star is slightly metal- 
rich compared to the Sun, in contrast to the photometric es- 
timate of [Fe/H] = -0.255 dex from the KIC. The initial set 
of spectroscopic constraints provided to the modeling teams 
(see ® came from a mean of the preliminary results from the 
analyses discussed above, with uncertainties large enough to 
cover the full range for each parameter: T e ff = 5635 ± 185 K, 
logg = 3.95 ± 0.25 dex, and [Fe/H] = -0.06 ± 0.25 dex. 

For the final estimate of the atmospheric parameters (see 
© we adopted the results from the VWA method, since it 
has been carefully tested again st direct methods fo r 10 nearby 
solar-type stars. Specifically, Brunt t et al.l (120 10a) used VWA 
to determine T e s from high-quality spectra and these values 
were compared to a direct method (nearly independent of 
model atmospheres) using the measured angular diameters 
from interferometry and bolometric flux measurements. A 
comparison of the direct (interferometric) and indirect (VWA) 
methods showed only a slight offset of -40 ± 20 K, and 
this offset h a s been removed for KIC 11026764. Similarly, 
iBruntt et"afl (1201 0a) determined the spectroscopic logg pa- 
rameter, which agrees very well for the three binary stars in 
their sample where the absolute masses and radii (and hence 
logg) have been measured. Based on these com parisons of 
direct and indirect methods. IBruntt et al.l (l2010al) also discuss 
the issue of realistic uncertainties on spectroscopic parame- 
ters, and their estimates of the systematic uncertainties have 
been incorporated into the adopted values from VWA listed in 
Table |2 



4. STELLAR MODEL SEARCH 

Starting with the minimal frequency set described in §|2] 
and the initial set of spectroscopic constraints from $3] eleven 
teams of modelers performed a "meta-search" of the parame- 
ter space. Each modeler had complete freedom to decide on 
the input physics and fitting strategy to optimize the match to 
the observations. The results of the individual fits were eval- 
uated in a uniform manner and ranked according to the total 
X 2 between the observed and calculated values of the indi- 
vidual oscillation frequencies and the spectroscopic proper- 
ties. These individual fits are listed in Table [3] and the details 
of the codes and fitting strategies employed by each team of 
modelers are described in the following subsections. 

4.1. Model A 

We employed the Aarhus st ellar evolution code (ASTEC; 
Christensen-Dalsgaard 2008a) for stellar evolution com- 
putations, and the ad i abatic pulsation package (ADIPLS; 
I Christensen-Dalsgaard! I2008bl) for frequency calculations. 
The input physics for the evoluti on calculations inclu ded the 
OPAL 2005 equation of state (iRoeers & Navfonovl 12002b . 
OPAL opacity tables (Iglesias & Rogers 1996) with low- 
temperature opacities from [Ferguson et al. (2005), and the 
NACRE nuclear reaction rates dAngulo et all 119991) . Con- 
vection was treated a ccording to the mixing-length theory of 
iBohm-Vitensd d 1958b . We did not include diffusion or con- 
vective overshoot in the models. 

We computed several grids of evolutionary tr acks spanning 
the pa rameter space around the values given by Chap lin et alJ 
(2010). We primarily adjusted the stellar mass and metal- 
licity in our grids, while fixing the mixing-length parameter 
a to 1.7. We scanned the parameter space in mass M from 
1.00 to 1.35 Mq\ initial heavy-element mass fraction Z; from 
0.009 to 0.025; and initial hydrogen mass fraction X\ from 
0.68 to 0.76. These values of Z\ and X\ cover a range of 
(Z/X)i = 0.012-0.037, or [Fe/H] = -0.317 to +0.177 dex us- 
ing [Fe/H] = log(Z/X)-log(Z/X) Q , wher e (Z/X) is the ratio 
at the stellar surface and (Z/X) Q = 0.0245 dGrevesse & Noels! 
1993). This range of [Fe/H] is compatible with the initial 
spectroscopic constraints. However, we later extended our 
grids to determine whether there was a better model with 
lower or higher metallicity. For all of the models on our 
tracks, we calculated the oscillation frequencies when the val- 
uesof r e ff and logg were within 2a of the derived values (see 
{33.31 l. We then assigned a goodness of fit to the frequency set 
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Table 3 

Initial model-fitting search for KIC 1 1026764. 



Model 


M/M 


Zs 


Y s 


a 


t (Gyr) 


L/Lq 


R/Rq 


7k(K) 




[Fe/H] 


x 2 


A.... 


1.13 


0.019 


0.291 


1.70 


5.967 


3.523 


2.026 


5562 


3.878 


-K).051 


9.8 


B.... 


1.23 


0.015 


0.250 


1.80 


4.861 


4.730 


2.097 


5885 


3.885 


-0.068 


12.4 


C... 


1.24 


0.021 


0.275 


1.79 


5.231 


4.293 


2.093 


5750 


3.890 


+0.120 


22.5 


D 


1.31 


0.044 


0.241 


1.42 


7.775 


2.598 


2.139 


5016 


3.895 


+0.400 


49.7 


E . . . . 


1.22 


0.013 


0.233 


1.83 


4.745 


4.710 


2.080 


5899 


3.888 


-0.100 


55.8 


F.... 


1.10 


0.014 


0.269 


1.85 


5.839 


3.850 


2.015 


5706 


3.870 


-0.102 


60.6 


G.... 


1.20 


0.022 


0.291 


1.88 


5.100 


3.960 


2.160 


5663 


3.880 


+0.116 


66.8 


H. ... 


1.10 


0.010 


0.250 


1.75 


6.752 


3.776 


2.012 


5677 


3.872 


-0.250 


284.3 


I 


1.27 


0.021 


0.280 


0.50 


3.206 


3.202 


1.792 


4854 


4.034 


+0.080 


637.3 


J 


1.12 


0.016 


0.276 


1.90 


6.505 


3.644 


2.026 


5593 


3.870 


-0.139 




zterror 


0.14 


0.007 


0.030 


0.78 


1.961 


0.678 


0.091 


188 


0.023 


±0.166 




K.... 


1.13 


0.017 


0.283 


1.80 


6.450 


3.610 


1.988 


5634 


3.890 


-0.044 




zterror 


0.13 


0.009 


0.009 




1.930 


0.770 


0.080 


161 


0.018 


±0.250 





of each model by calculating x 2 : 



n.l 



v°^(n)-v™ oae \n) 
cr(V° bs (n)) 



(2) 



where N is the total number of modes included, v° bs (n) and 
vf° Ael (n) are the observed and model frequencies for a given 
spherical degree I and radial order n, while a{vf b!i {n)) repre- 
sents the uncertainties on the observed frequencies. We cal- 
culated Y 2 after correcting the f requencies for surface effects, 
following lKjeldsen et al.l d2008h : 



^ob S («)-^best(«) = a 



*W«) 
MO 



(3) 



where v <a s {n) and twin) are the observed and best model 
frequencies with spherical degree I = and radial order n, 
and isq is the frequency of maximum power in the oscilla- 
tion spectrum, which is 857/iHz for KIC 1 1026764. We fixed 
the exponent b to the value derived for the Sun (b = 4.90) by 
iKjeldsen et all (120081) . and a was calculated for each model. 
We computed smaller and more finely sampled grids around 
the models with the lowest \ 2 to refine the fit. The prop- 
erties of the best model are listed in Table [3] Although we 
found models with higher or lower logg that had large sepa- 
rations quite close to the observed value, the individual fre- 
quencies were not close to the observations, and the larger 
range of metallicity did not yield improved results. We also 
found more massive models (around 1.3 M Q ) with a total x 2 
value comparable to our best fit, but they did not include the 
mixed modes. Model A is the best match from the family of 
solutions (also including Models F, H, J and K) with masses 
near 1.1 M®. 

4.2. Model B 

W e used the Geneva stellar evolution code including rota- 
tion {Eggenberger et al. 2008) for all co mputations. This code 
includ es the OPAL equatio n of state dRogers & Navfonovl 
120021) , the OPAL opacities dlglesias & Rogers! 1 19961 ) com- 
ple mented at low tempera t ures w ith the molecular opacities 
of lAlexander & Ferguson! d 19941) . the NACRE nuclear re- 
letalJ I 



action rates dAngulo i 



and the standard mixing- 



length formalism for convection dBohm- Vitensd 19 58). Over- 
shooting from the convective core into the surrounding ra- 
diatively stable layers by a distance d m = a ov mm[H p ,r coie ] 
(Maeder & Mevnet 1989) is included with an overshoot pa- 
rameter a ov = 0.1. 

In the Geneva code, rotational effects are computed in the 
framework of shellular rotation. The transport of angula r mo- 
mentum then obeys an adv ection-diffusion equation dZahnl 
1 19921 iMaeder & Zahnlll998l) . while the vertical transport of 
chemicals through the combined action of vertical advection 
and strong horizo ntal diffusion can be des cribed as a purely 
diffusive process dChabover & Zahijl 1992b . Since the model- 
ing of these rotational effects has be en described in previous 
papers (e.g. Eggenberger et al.ll2010l) . we simply note that the 
Geneva code includes a comprehensive treatment of shellular 
rotation and that meridional circulation is treated as a truly 
advective process. For a detailed analysis of the effect of cen- 
trifugal force on the oscillation frequencies, see Appendix lAl 
In addition to rotation, atomic diffusion of helium and heavy 
elements is included with diffu sion coefficients cal culated ac- 
cording to the prescription of Paq uette et al.l d!986l) . 

The properties of a stellar model including rotation depend 
on six parameters: the mass M, the age t , the mixing-length 
parameter a = l/H p for convection, the initial rotation veloc- 
ity on the ZAMS and two parameters describing the initial 
chemical composition of the star. For these two parameters, 
we chose the initial helium abundance Yj and the initial ratio 
between the mass fraction of heavy elements and hydrogen 
{Z/X)\. This ratio can be related to the metallicity [Fe/H] 
assuming that \og{Z/X) = [Fe/H] + log(Z /X) w ; we adopt 
the so lar value (Z/X)q = 0.0245 given by iGrevesse & Noelsl 
( 1993). For these computations, the mixing-length parameter 
was fixed to a solar calibrated value (q = 1 .7998) and the ini- 
tial rotatio n velocity on th e ZAMS was 50kms _1 . The brak- 
ing lawofKawiiiJilli) was used to reproduce the magnetic 
braking experienced by low-mass stars during main-sequence 
evolution. 

With the above assumptions, the characteristics of a stel- 
lar model depend on only four parameters: M, f, Y\ and 
(Z/X)i. The determination of the parameters that best repro- 
duce the observationa l constraints was then pe r forme d in two 
steps as described in Eggenberge r & Carried d2006l) . First, 
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a grid of models with global properties in reasonable agree- 
ment (within 2a) with the adopted spectroscopic constraints 
was constructed. Theoretical frequencies of / < 2 modes in 
the observed range of 590-1 100 mHz were computed using 
the adiabatic pulsation code dChristensen-D alsgaard 2008b) 
along with the characteristic frequency separations. The mean 
large separation was determined by considering only radial 
modes. The effects of incomplete modeling of the external 
layers on computed frequencies were taken into account us - 
ing the empirical power law given by Kjeld sen et al.1 (2008). 
This correction was applied to theoretical frequencies using 
the solar calibrated value of the exponent (b = 4.90) and cal- 
culating the coefficient a for each stellar model. 

Using spectroscopic measurements of [Fe/H], r e ff and logg 
together with the observed frequencies, a x 2 minimization 
was performed to determine the set of model parameters that 
resulted in the best agreement with all observational con- 
straints. The properties of the best model are listed in Ta- 
ble[3] This model correctly reproduces the spectroscopic mea- 
surements of the surface metallicity and logg, but exhibits a 
slightly higher effective temperature. It is in good agreement 
with the asteroseismic data and in particular with the observed 
deviation of the I = 1 modes from asymptotic behavior. Model 
B is the best match from the family of solutions (also includ- 
ing Models C, E and G) with masses near 1.2 M Q . 

4.3. Model C 

The Garching Ste llar Evolution Code (GARSTEC; 
iWeiss & Schlattlll2008l) is a one-dimensional hydrostatic code 
which does not include the effects of rotation. For the model 
calcu lations we used the OPAL equation of state (Roge rset al.l 
1996) complem ented with the MHD equ ation of state at low 
temperatures dHummer & Mihalas 1988), OPAL opacities for 
high temperatures (Iglesias & Rogers 1996) and Ferg uson's 
opacities for low t emper atures dFerguson et al.l 120051) . the 
iGrevesse & Sauvall (119981) solar mixture, and the NACRE 
comp ilation of thermonuclear reaction rates d An gulp et alj 
1999). Mixing is performed diffusively in convective regions 
using the m ixing-length theory for convect ion in the formu- 
lation from Kippenhahn & Weigert] (119901) . and convective 
overshooting can optionally be implemented as a diffusive 
process with an exponential decay of the convective velocities 
in the radiative zone. The amount of mixing for overshooting 
depends on an efficiency parameter A calibrated with open 
clusters (typically A = 0.016). A tomic diffusion ca n be 
applied following the prescription of T houl et alJ (119941) . and 
we use a plane-parallel Eddington grey atmosphere. 

We started all of our calculations from the pre-main se- 
quence phase. The value of the mixing-length parameter 
for convection was fixed (a = 1.791 from our solar calibra- 
tion), the Schwarzschild criterion for definition of convec- 
tive boundaries was used, and we did not consider convec- 
tive overshooting or atomic diffusion. We constructed a grid 
of models in the mass range between 1 .0 M Q and 1 .3 M Q (in 
steps of 0.01) for several [Fe/H] values from the spectroscopic 
analysis: 0.06, 0.09, 0.12, and 0.15. To convert the observed 
values into total metallicity, we applied a chemical enrichment 
law of AY I AZ = 2 and used the primordial abundances from 
our solar calibration. We did not explore variations in either 
the hydrogen or helium abundances. 

Once all of the tracks were computed, we restricted our 
analysis to those models contained within the spectroscopic 
uncertainties. For these cases, we calculated the oscillation 
frequencies using the adiabatic pulsation package (ADIPLS; 



IChristensen-Dalsgaardl 12008b!) and looked for the model 
which best reproduced the large frequency separation (no sur- 
face correction was applied to the calculated frequencies). 
Several models were found to fulfill these requirements for 
each metallicity grid, and among those best fit models to the 
large frequency separation we then performed a \ 2 test to ob- 
tain the global best fit to the individual frequencies and the 
spectroscopic constraints. Our global best fit model came 
from the grid with [Fe/H] = 0.12, and the properties are listed 
in Table [3] This model reproduces well the observed mixed 
modes. 

4.4. Model D 

The Asteroseismic Modeling Portal (AMP) is a 
web-based tool tied to TeraGrid computing resources 
that uses the Aarhus stell ar evolution code (ASTEC; 
Christensen-Dalsgaard 2008a) and adiabatic pulsation code 
(ADIPLS; [Christensen-Dalsgaard 2008b) in conjunction 
with a parallel genetic algorithm (Metcalfe & Charb onneaul 
1 20031) to optimize the match to observational data (see 
Metcalfe et ail 120091) . The models use the OP AL 2005 
equation of state (see [Rogers & Navfonovs 2002) a nd th e 
most recent OPAL opacities (see Iglesias & Ro gers! 119961) . 
supplemented by Kurucz opacities a t low temperatures. The 
nuclea r reaction rates come from iBahcall & Pinsonneaultl 
d 19951). convection is de scribed by the mixing-length theory 
of iBohm-Vitensel £1958), and we can op tionally include the 
effects of helium settling as described by Michau d & Proffittl 
(fl99l . 

Each model evaluation involves the computation of a stel- 
lar evolution track from the zero-age main sequence (ZAMS) 
through a mass-dependent number of internal time steps, ter- 
minating prior to the beginning of the red giant stage. Rather 
than calculate the pulsation frequencies for each of the 200- 
300 models along the track, we exploit the fact that the aver- 
age frequency separation of consecutive radial orders (Ai/q) 
in most cases is a monotonic ally decreasing function of age 
dChristensen-Dalsgaardl [T993 ). Once the evolution track is 
complete, we start with a pulsation analysis of the model at 
the middle time step and then use a binary decision tree — 
comparing the observed and calculated values of (Av^) — to 
select older or younger models along the track. This allows 
us to interpolate the age between the two nearest time steps 
by running the pulsation code on just 8 models from each 
stellar evolution track. The frequencies of each model are 
th en corrected for surfa ce effects following the prescription 
of lKjeldsen eTai] d2008l) . 

The genetic algorithm (GA) optimizes four adjustable 
model parameters, including the stellar mass (M) from 0.75 
to 1.75 M Q , the metallicity (Z) from 0.002 to 0.05 (equally 
spaced in logZ), the initial helium mass fraction (Yj) from 0.22 
to 0.32, and the mixing-length parameter (a) from 1 to 3. The 
stellar age (f ) is optimized internally during each model eval- 
uation by matching the observed value of (Ai/q) (see above). 
The GA uses two-digit decimal encoding, such that there are 
100 possible values for each parameter within the specified 
ranges. Each run of the GA evolves a population of 128 mod- 
els through 200 generations to find the optimal set of parame- 
ters, and we execute 4 independent runs with different random 
initialization to ensure that the best model identified is truly 
the global solution. The resulting properties of the optimal 
model are listed in Table [3] 

The extreme values in this global fit arose from treating 
each spectroscopic constraint as equivalent to a single fre- 
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quency. Since the adopted spectroscopic errors are large, they 
provide much more flexibility for the models compared to the 
individual frequencies with relatively small errors. Conse- 
quently, the fitting algorithm found it advantageous to shift the 
effective temperature and metallicity by several a from their 
target values to achieve significantly better agreement with the 
22 oscillation frequencies. The improvement in the frequency 
match outweighed the degradation in the spectroscopic fit for 
the calculation of \ 2 - The solution to this problem may be 
to calculate a separate value of \ 2 f° r the asteroseismic and 
spectroscopic constraints, and then average them to provide 
more equal weight to the two types of constraints. This is an 
important lesson for future automated searches, and explains 
why Model D does not align with either of the two major fam- 
ilies of solutions. 

4.5. Model E 

We used the Yal e Stellar Evolution Code (YREC; 
iDemarque et al.l 120081) in its non-rotating configuration to 
model KIC 11026764. All models were constructe d with 
the OPAL equation of state dRogers & NaYfonovll2002l). We 
used OPAL high temperature opacities dlglesias& Rogers 
19961) supplemented with low temperature opacities from 
Ferguson et al. (2005). The NACRE nuclear reaction rates 
( Angu kTet al.l l 1999) were used. We assumed that the curren t 
solar metallicity is that given by iGrevesse & Sauvall (119981) . 
We have not explored the c onsequences of using the lower 
metallicity measurements of Aspl und et al.1 (120051) or the in - 
termediate metallicity measurements of iLudwig et al.l d2009t) . 
We searched for the best fit within a fixed grid of models. 
There were eight separate grids defined by different combina- 
tions of the mixing-length parameter (a = 1.83 or a = 2.14), 
the initial helium abundance (either F; = 0.27 or Yj calculated 
assuming a AY / AZ = 2, with Yj for [Fe/H] = being the cur- 
rent solar CZ helium abundance), and the amount of over- 
shoot (0H p or Q.2H p ). All models included gravitational set- 
tling of helium an d heavy elements using the formulation of 
iThoul et all (fl99l . 

Our fitting method included two steps. In the first step, 
we calculated the average large frequency separation for the 
models and selected all of those that fit the observed sepa- 
ration within 3cr errors. We adopted the obser ved value of 
Av = 50.8 ± 0.3^Hz from IChaplin et all (120101) . A second 
cut was made using the effective temperature: all models 
within ±2ct of the observed value were chosen. A third cut 
was made using the frequencies of the three lowest-frequency 
1 = modes. Given the small variation in mass, this process 
was effectively a radius cut. The selected models had radii 
around 2R Q . Note that the Yale-Birmingham radius pipeline 
dBasu etalJ2010l) finds a radius of 2. 1 8!g jQR Q for this star us- 
ing the adopted values of Av, T e fi, logg and [Fe/H]. In the sec- 
ond step of the process, we made a finer grid in mass and age 
around the selected values and then compared the models with 
the observed set of frequencies. The properties of our best fit 
model are listed in Table [3] This model was constructed with 
Y, = 0.27, Z = 0.0147 and core overshoot of 0.2H p . We were 
unable to find a good model without core overshoot. 

4.6. Model F 

We modeled KIC 1 1026764 with the Catani a Astrophysical 
Obser vatory version of the GARSTEC code (Bonan no et alj 
2002) using a grid-based approach. The input physics of 
this stellar evolution code included the OPAL 2005 equa- 



tion of state dRogers & Navfonovll2 002) and the OPAL opac- 
ities (Igl esias & Rogersll 19961) comple mented in the low tem- 
peratu re regime with the tables of lAlexander & Fer guson 
(1994). The nuclear r eaction rates were taken from the 
NACRE collaboration (Anguloetal. 1999) and the stan- 
dard mixing-length formalis m for convection was used 
(Kippenhahn & Weigert 1990). Microscopic diffusion of hy- 
drogen, helium and all of the major metals can optionally be 
taken into account. The outer boundary conditions were de- 
termined by assuming an Eddington grey atmosphere. 

A grid of evolutionary models was computed to span the la 
uncertainties in the spectroscopic constraints obtained from 
ground based observations. When a given evolutionary track 
was in the error box, a maximum time step of 20 Myr was 
chosen and the frequencies were computed with the ADIPLS 
code. A non-uniform grid of mass in the range l.O-1.24M , 
helium abundance in the range Y\ = 0.26-0.31, mixing-length 
parameter a = 1.6-1.9 and initial surface heavy-element abun- 
dances (Z/X)i = 0.022-0.029 was scanned. A global optimiza- 
tion strategy was implemented by minimizing the x 2 for all of 
the I = 0, / = 1 and 1 = 2 modes . The empirical surface effect, 
as discussed by Kiel dsen et al.l d2008l) . was used to correct all 
theoretical frequencies. The properties of the best model with 
heavy-element diffusion and surface corrected frequencies are 
listed in Table [3] 

4.7. Model G 

We use d a version of the Aarhus st ellar evolution code 
(ASTEC; Chris tensen-Da lsgaard 2008a ) which incl udes the 
OPAL 20 01 equation of state dRogers et all 119961). OPAL 
opacit ies dlglesias & Rogers] 119961) . iBahcall & PinsonneaujJ 
d!995|) nuclear cross sect ions and the mixing-length formal- 
ism dBohm-Vitensd[l958l) for convection. We computed sev- 
eral grids of models by varying all of the input parameter s 
within the range of the observed errors (Chaplin etal. 2010). 
In particular we calculated evolutionary tracks by varying 
the input mass in the range M = 0.9-1.2 M Q , the metallic- 
ity in the range Z = 0.009-0.03, and the hydrogen abundance 
in the range X = 0.67-0.7. We also adopted different val- 
ues of the mixing-length parameter in the range a = 1.67- 
1.8 8. We calculated addition al evolutionary models using 
the ICanuto & Mazzitelli ( 1992) convection formulation. The 
mixing-length parameter a of the CM formulation was chosen 
in the range a = acM = 0.9-1.0. To obtain the deviation from 
asymptotic behavior observed in the / = 1 modes (the mixed 
modes) we did not include overshooting in t he cal culation, 
following the conclusion of iDi Mauro et al.l (120031) . These 
models are distinct from the grid used to produce Model A, 
not only because they employ a slightly older EOS and nu- 
clear reaction rates, but also because the grid search included 
a and calculated fewer models within the specified range 
of paramet er values. We used the adiabatic oscillation code 
(ADIPLS; Christensen-Dalsg aardll2008bl) to calculate the p- 
mode eigenfrequencies with harmonic degree / = 0-2. The 
characteristics of the model which best fits the observations 
are listed in Table [3] 

4.8. Model H 

For the evolution calculations we used the publicly avail- 
able Dartmouth stellar evolut i on code (DSEP;IChabover et all 
l200lt iGuenther et all W99% iDotter et all |2007|), which~~is 
based on t he code developed by Pie r re Demarque and his 
students (lLarson & D emarquel 119641; [Demarque & Mengel 
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1197 ll) . The i nput physics includes h igh temperature opacities 
from OPAL (Iglesia s & Rogers] | 19961) low temperature opac- 
iti es from [F erguson et al.l (l2005h . the nuclear reaction rates 
of Bahcall & Pi nsonneaultl (1199 2). helium and h eavy-element 
settling and diffusion (Michau d & Proffittll 19931). and Debye- 
Hiick el corrections to the equation of state (Guenther et al. 

1992) . The models employ the standard mixing-length the- 
ory. Convective core overshoot is calculated assuming that 
the extent is proportional to the pressure scale height at 
the boundary (Demarque et al. 2004). The models used the 
standard conversion fr om [Fe/H] and AF/AZ to Z and Y 
(Chab oyer et al.l fl999l) . The oscillation frequ encies were 
compu ted using the adiabatic oscillation co des of lKosovichevI 
(1999) and lChristensen-Dalsgaardl (l2008bl) . No surface cor- 
rections were applied. 

The strategy to find a model matching the observed spec- 
troscopic constraints involved calculating a series of evolu- 
tionary tracks in the logg-T^ff plane for a mass range of 1.0- 
1.3 Mq, heavy-element abundance Z = 0.01-0.03, initial he- 
lium abundance Y\ = 0.25-0.30, and mixing-length parameter 
a = 1.70-1.75. We then selected the models closest to the tar- 
get values within half of the specified uncertainties. All mod- 
els for the search were calculated assuming an overshoot pa- 
rameter a ov = 0.2, and included element diffusion. For com- 
parison, the corresponding models without diffusion and con- 
vective overshoot were also calculated. The oscillation spec- 
tra were matched to the observed frequencies, first by com- 
paring the frequencies of radial (/ = 0) modes with the cor- 
responding observed frequencies, and then selecting a model 
with the closest frequency values for the / = 1 and / = 2 modes. 
The properties of the final model are listed in Table [3] This 
model matches the observed frequencies quite well except for 
the first two I = 1 modes, which deviate by about 12-13 /iHz. 
Our search demonstrated that the behavior of the mixed mode 
frequencies is sensitive to element diffusion and convective 
overshoot. This requires further investigation. 

4.9. Model I 

To characterize KIC 11026764 we const ructed a gri d 
of stellar models with the CESAM code (iMorell 11997b . 
and computed their oscillati on frequencies with the adia- 
batic oscillation code FILOU ( Suaiez l2002l:LSuarez&Goupill 
12008b . Opacity tables were taken from the OPAL 
package (Iglesias & Rogers 1996), complemented at low 
temperatures (T < 1Q 4 K) w ith the tables provided by 
lAlexandeT & Ferguson (1994). The atmosphere was con- 
structed from a Eddington T-t relation and was assumed to 
be grey. The stellar metallicity (Z/X) wa s derived from the 
JTe/H] value assuming (Z/X) Q = 0.0245 (IGrevesse & Noels! 

1993) , T pl = 0.235 and Z pr = for the primordial helium and 
heavy-element abundances, and a value AF/AZ = 2 for the 
enrichment ratio. No microscopic diffusion of elements was 
considered. 

The main strategy was to search for representative equilib- 
rium models of the star in a database of 5 x 10 5 equilibrium 
models, querying for those matching the global properties 
of the star, including the effective temperature, gravity and 
metallicity. Using this set of models, we then applied the as- 
teroseismic constraints, including the individual frequencies 
and large separations. The global fitting method involved a 
X 2 minimization, taking into account all of the observational 
constraints simultaneously. No correction for surface effects 
was applied, and no a priori information on mode identifica- 



tion was assumed when fitting the individual frequencies. The 
properties of the best model we found is listed in Table[3] and 
includes overshooting with a ov = 0.3. Note that the analysis 
did not adopt the identifications of spherical degree (/) from 
^2] and it included I = 3 modes to perform the match. As a 
consequence, the final result is much different than any of the 
others and it does not fall into either of the two major families 
of solutions. 

4.10. Model] 

The SEEK procedure makes use of a large grid of 
stellar models comp uted with the Aarhus stellar evolu- 
tion code (ASTEC; IChristensen-Dalsgaardl l2008al) . It 
compares the observations with every model in the grid 
and makes a probabilistic assessment, with the help of 
Bayesian statistics, about the global properties of the star. 
The model grid includes 7,300 evolution tracks contain- 
ing 5,842,619 individual models. Each track begins at the 
ZAMS and continues to the red giant branch or a max- 
imum age of t = 15 Gyr. The tracks are separated into 
100 subsets with different combinations of metallicity Z, 
initial hydrogen content X\ and mixing-length parameter 
a. These combinations are separated into two regularly 
spaced and interlaced subgrids. The first subgrid comprises 
tracks with Z = [0.005,0.01,0.015,0.02,0.025,0.03], X { = 
[0.68,0.70,0.72,0.74], and a = [0.8,1.8,2.8] while the sec- 
ond subset has Z = [0.0075,0.0125,0.0175,0.0225,0.0275], 
Xi = [0.69,0.71,0.73], a = [1.3,2.3]. Every subset is com- 
posed of 73 tracks with masses between 0.6 and 3.0 Mq. The 
spacing in mass between the tracks is 0.02 Mq from 0.6 to 
1.8 Mq and 0.1 from 1.8 to 3.0 M . A relatively high value 
of Yq = 0.2713 and Z© = 0.0196 for the Sun has been used 
for the standard definition of [Fe/H] in SEEK. This value is 
used to cali brate solar models from AS TEC to the correct lu- 
minosity (Chris tensen-Dalsgaardll9 98). The input physics in- 
clude the OPAL equation of s tate, opacity tables from OPAL 
(Iglesias & Rogers 1996) a nd [Alexander & Fergusonl (fl994). 
and the metallic mixture of IGrevesse & Sauv al ( 1998). Con- 
vection is trea t ed ac cording to the mixing-length theory of 
iBohm-Vitensd dl958l) with the convective efficiency charac- 
terized by the mixing-length to pressure height scale ratio a, 
which varies across the grid of models. Diffusion and over- 
shooting were not included. 

The grid allows us to map the physical input parameters 
of the model p = {M,t,Z,X[,a} into the grid of observable 
quantities q g = {A^,<5^,r e ff,logg, [Fe/H],...}, defining the 
transformation 

q g = /C(p). (4) 

We compare these quantities to the observed values q obs with 
the help of a likelihood function £, 

£ -(nds;) (s 

and the usual \ 2 definition 




where er, is the estimated error for each observation q° , and 
N is the number of observables. The maximum likelihood is 
then combined with the prior probability of the grid /o to yield 
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Table 4 

Final model-fitting results for KIC 1 1026764. 



Model 


M/Mp, 


z s 


Y s 


a 


t (Gyr) 




R/R G 


7k(K) 


logg 


[Fe/H] 


x 2 


FA... 


1.13 


0.017 


0.305 


1.64 


5.268 


4.141 


2.036 


5778 


3.872 


+0.009 


3.69 


AA .. 


1.13 


0.019 


0.291 


1.70 


5.977 


3.520 


2.029 


5556 


3.877 


+0.051 


6.11 


AA\. 


1.13 


0.019 


0.291 


1.70 


5.935 


3.454 


2.026 


5534 


3.877 


+0.031 


7.40 


GA .. 


1.10 


0.017 


0.296 


1.88 


6.100 


3.420 


2.010 


5539 


3.870 


+0.004 


78.05 


CA... 


1.13 


0.019 


0.291 


1.70 


6.204 


3.493 


2.030 


5546 


3.876 


+0.050 


152.91 


EA... 


1.12 


0.019 


0.291 


1.70 


6.683 


3.202 


2.029 


5424 


3.870 


+0.050 


230.58 


AB. 


1.23 


0.018 


0.242 


1.80 


5.869 


3.804 


2.083 


5591 


3.890 


-0.010 


6.97 


AB\. 


1.20 


0.024 


0.276 


1.80 


5.994 


3.460 


2.072 


5475 


3.884 


+0.146 


7.26 


BB... 


1.22 


0.021 


0.270 


1.80 


5.153 


4.190 


2.061 


5758 


3.896 


+0.072 


7.57 


FB... 


1.24 


0.021 


0.280 


1.79 


4.993 


4.438 


2.092 


5800 


3.890 


+0.091 


8.52 


EB... 


1.22 


0.013 


0.232 


1.80 


4.785 


4.651 


2.079 


5882 


3.890 


-0.130 


18.54 


CB... 


1.24 


0.015 


0.250 


1.80 


5.064 


4.696 


2.089 


5887 


3.892 


-0.080 


45.84 


J' .... 


1.27 


0.021 


0.270 


1.52 


4.260 


4.011 


2.105 


5634 


3.892 


+0.080 




± error 


0.09 


0.003 


0.024 


0.74 


1.220 


0.371 


0.064 


81 


0.020 


±0.060 




K' ... 


1.20 


0.022 


0.278 


1.80 


5.980 


3.700 


2.026 


5619 


3.900 


+0.070 




± error 


0.04 


0.003 


0.003 




0.610 


0.300 


0.027 


79 


0.006 


±0.060 





the posterior, or the resulting probability density 

/(p)cx/o(p)£(/C(p)). (7) 

This probability density can be integrated to obtain the value 
and uncertainty for each of the parameters, as listed in Ta- 
bles|3]and|4] It can also be projected onto any plane to get the 
correlation between two parameters, as shown in AppendixlBl 
The details of the SEEK procedure, including the choice of 
pr iors, and an introduction to Bayesian statistics can b e found 
in Oui rion. Christensen-Dalsgaard & Arentoftl (120101) . 

SEEK uses the large and small separations and the median 
frequency of the observed modes as asteroseismic inputs. Val- 
ues of A^o = 50.68 ± 1 .30 (computed with / = modes only) 
and (5^o.2 = 4.28 ± 0.73 (computed from I = 0, 2 modes) were 
derived from Kepler data around a central value 900 /iHz. 
Each separation is the mean of the individual observed separa- 
tions, while the error is the standard deviation of the individ- 
ual values from the mean. These v alues differ slightly from 
those given in IChaplin et alj (2010) because they are calcu- 
lated from the individual frequencies rather than derived from 
the power spectrum. Using these asteroseismic inputs along 
with the initial spectroscopic constraints, we obtained the pa- 
rameters listed in Table [3] For a SEEK analysis of the impor- 
tance of the asteroseismic constraints, see AppendixlBl 

4.11. Model K 

To investigate how well we can find an appropriate model 
without comparing individu al oscillation fre quencies, we 
used the RADIUS pipeline (IStello et al.ll2009al) . which takes 
T e ff, \ogg, [Fe/H], and Av as the only inputs to find the 
best fitting m o del. T he value of Av = 50.8 ± 0.3 fiHz from 
ChaplinetalJ d20Toh was adopted. The pipeline is based 
on a l arge grid of ASTEC models (IChristensen-Dal sgaard 
l2008al) using the EFF equation of state ([E ggleton et al.ll 19731) 7 
We use the opacity tables of Ro gers & Iglesiasl (11995b and 
Kuruczl d 199 11) for T < 10 4 K with the solar mixture of 
Grevesse & Noelsl (119931) . Rotation, overshooting and diffu- 
sion were not included. The grid was created with fixed values 
of the mixing-length parameter (a = 1.8) and the initial hy- 



drogen abundance (Xi = 0.7). The resolution in logZ was 0.1 
dex between 0.001 < Z < 0.055, and the resolution in mass 
was 0.01 M Q from 0.5 to 4.0 M Q . The evolution begins at 
the ZAMS and continues to the tip of the red giant branch. 
To convert between the mode l values of Z and the observed 
[Fe/H], we used Z Q = 0.0188 (ICoxll2000l) . 

We made sli ght modifications t o the RADIUS approach 
described by IStello et all (l2009al) . First, the large fre- 
quency separation was der ived by scaling the solar value 
( Kield sen & Bed ding 1995) instead of calculating it directly 
from the model frequencies. Although there is a known 
systematic difference between these two ways of deriving 
Ay, the effect is p robably below the 1% level (Ste llo et al.l 
2009b; iBasu et alj |2010). Second, we pinpointed a single 
best-fitting model based on a \ 2 formalism that was applied 
to all models within ±3<7 of the observed properties. The 
properties of the best fitting model are listed in Table [3] This 
model shows a frequency pattern in the echelle diagram that 
looks very similar to the observations if we allow a small 
tweaking of the adopted frequency separation (52.1 /iHz) 
used to generate the echelle. This basically means that we 
found a model that homologously represents the observations 
quite well. In particular, we see relative positions of the 
/ = 1 mode frequencies that are very similar to those observed. 



5. MODEL-FITTING RESULTS 

Based on the value of x 2 from the initial search in $4] we 
adopted two reference models (Models A and B in Table [3j> 
each constructed with a very different set of input physics, 
but almost equally capable of providing a good match to the 
observations. Note that Model A does not include overshoot 
or diffusion, while Model B includes overshoot, diffusion and 
a full treatment of rotation. These two models differ signif- 
icantly in the optimal values of the mass, effective temper- 
ature, metallicity and luminosity, but they both agree with 
the observational constraints at approximately the same level 
(\ 2 ~ 10). With the exceptions of Models D and I (see sub- 
sections above), the other independent analyses generally fall 
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into the two broad families of solutions defined by Models A 
and B. The lower mass family includes Models A, F, H, J and 
K, while the higher mass family includes Models B, C, E and 
G. We identified several additional asteroseismic constraints 
from the maximal frequency set (see © and we adopted re- 
vised spectroscopic constraints from VWA (see ® to refine 
our analysis of Models A and B using several different codes. 

5.1. Refining the Best Models 

Comparing the theoretical frequencies of Models A and B 
with the maximal frequency set from $2] we identified four 
of the seven additional oscillation modes that could be used 
for refined model-fitting (see Table [TJ. Recall that the maxi- 
mal frequency set comes from the individual analysis with the 
smallest rms deviation with respect to the maximal list, so the 
frequencies of the modes from the minimal set are slightly 
different in the maximal set. Without any additional fitting, 
these subtle frequency differences improve the \ 2 of Models 
A and B when comparing them to those modes from the max- 
imal set that are also present in the minimal set. There is one 
additional I = mode (n = 12) and three additional / = 2 modes 
(n = 11,16, 19) in the maximal set that are within 3er of fre- 
quencies in both Models A and B. Considering the very differ- 
ent input physics of these two models, we took this agreement 
as evidence of the reliability of these four additional frequen- 
cies and we incorporated them as constraints for our refined 
model-fitting. Two of the remaining frequencies in the maxi- 
mal set (n = 12, 1 = 1 and 2) were not present in either Models 
A or B, while one (n = 10, 1 = 2) had a close match in Model B 
but not in Model A. We excluded these three modes from the 
refined model-fitting. Given that the / = 1 mo des provide the 
strongest constraints on the models (see 35.2b . the additional 
I = and 1 = 2 modes are expected to perturb the final fit only 
slightly. 

In addition to the 26 oscillation frequencies from the maxi- 
mal set, we also included stronger spectroscopic constraints in 
the refined model-fitting by adopting the results of the VWA 
analysis instead of using the mea n atm ospheric parameters 
from the preliminary analyses (see 33.3b . Although the uncer- 
tainties on all three parameters are considerably smaller from 
the VWA analysis, the actual values only differ slightly from 
the initial spectroscopic constraints. These were just three of 
the 25 constraints used to calculate the x 2 an d ran k the ini- 
tial search results in Table [3] Since the 22 frequencies from 
the minimal set were orders of magnitude more precise, they 
dominated the x 2 determination. Although the spectroscopic 
constraints from VWA are more precise than the initial atmo- 
spheric parameters, they are still much less precise than the 
frequencies and should only perturb the x 2 ranking slightly. 
Consequently, we do not need to perform a new global search 
after adopting the additional and updated observational con- 
straints. 

Several modeling teams used the updated asteroseismic and 
spectroscopic constraints for refined model-fitting with a vari- 
ety of codes. Each team started with the parameters of Models 
A and B from Table [3] and then performed a local optimiza- 
tion to produce the best match to the observations within each 
family of solutions. The results of this analysis are shown in 
Table [4] where the refined Models A and B are ranked sep- 
arately by their final x 2 value. Each model is labeled with a 
letter from Table [3] to identify the modeling team, followed 
by either A or B to identify the family of solutions. The two 
pipeline approaches labeled J' and K' simply adopted the re- 
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Figure 3. An echelle diagram showing the 26 frequencies from the maxi- 
mal set that were used as constraints (solid points) with the frequencies of 
Models AA (blue) and AB (red) for comparison. The star exhibits modes 
with / = (circles), 1=1 (triangles), and 1 = 2 (squares). A greyscale map 
showing the power spectrum (smoothed to 1 /xHz resolution) is included in 
the background for reference. 

vised constraints to evaluate any shift in the optimal param- 
eter estimates and errors. Note that both the SEEK and RA- 
DIUS pipelines identified parameters in the high-mass family 
of solutions when using the revised constraints, but the low- 
mass family was only marginally suboptimal. The apparent 
bifurcation of results in Table H] into two values of a arises 
from the decision of most modelers to fix this parameter in 
each case to the original value from Table [3] For each fam- 
ily of solutions, the modeling teams adopted the appropriate 
input physics: neglecting overshoot and diffusion for the re- 
fined Models A, while including both ingredients for the re- 
fined Models B. One team produced two additional models 
(labeled AA' and AB') to isolate the effect of input physics 
on the final results. Model AA' started from the parameters 
of Model A but included overshoot and helium settling, while 
Model AB' searched in the region of Model B but neglected 
overshoot and diffusion. 

5.2. Stellar Properties & Error Analysis 

An inspection of the results in Table |4]reveals that models 
in either family of solutions can provide a comparable match 
to the observational constraints. This ambiguity cannot be at- 
tributed to the input physics, since the models that sample all 
four combinations of the input physics and family of solutions 
(AA, AA', AB, AB') have comparable x 2 values. The indi- 
vidual frequencies of these models all provide a good fit to 
the data, including the / = 1 mixed modes. The observations 
are compared to two representative models in Figure [3] us- 
ing an echelle diagram, where we divide the frequency spec- 
trum into segments of length ( Ais) and plot them against the 
oscillation frequency. This representation of the data aligns 
modes with the same spherical degree into roughly vertical 
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Figure 4. The significance S, of each observational constraint in determining the parameters of Model AA. Each oscillation frequency is labeled with the 
reference radial order from Table[T] and the degrees / = 0, 1 and 2 are shown sequentially in different shades from left to right. 



columns, with / = modes shown as circles, / = 1 modes 
shown as triangles, and / = 2 modes shown as squares. The 
26 modes from the maximal frequency set that were included 
in the final fit are indicated with solid points. Open points in- 
dicate the model frequencies, with Model AA shown in blue 
and Model AB shown in red. A greyscale map of the power 
spectrum (smoothed to 1 /iHz resolution) is included in the 
background for reference. The different models appear to be 
sampling comparable local minima in a correlated parameter 
space. Without additional constraints, we have no way of se- 
lecting one of these models over the other. 

We can understand the two families of models by con- 
sidering the general properties of subgiant stars, where a 
wide range of masses can have the same stellar luminos- 
ity with minor adjustments to the input physics and other 
model parameters — in particular the helium mass fraction. 
This degeneracy between mass and helium abundance has 
been discussed in the modelin g of specific subgiant stars (e.g . 
Fernandes & Monteiro 2003; Pinheiro & Fernandes 2010; 
Yang & Meng 2010), and it adds a large uncertainty to 



the already difficult problem of dete rmining the helium 
abundance in main-seque nce stars (e.g. iVauclair et all 120081 ; 
iSoriano & Vauc lair 2010). The luminosity of stellar models 
on the subgiant branch is mainly determined by the amount 
of energy produced at the edge of the helium core established 
during the main-sequence phase. The rate of energy produc- 
tion depends on the temperature and the hydrogen abundance 
in that layer. As a consequence it is possible to find models at 
the same luminosity with quite different values of total mass, 
by adjusting the other parameters to yield the required temper- 
ature at the edge of the helium core. Fortunately, for subgiant 
stars the presence of mixed modes can provide additional con- 
straints on the size of the helium core. For KIC 1 1026764 this 
mixed mode constraint significantly reduces the range of pos- 
sible masses to two specific intervals, where the combination 
of stellar mass and core size are compatible with the atmo- 
spheric parameters and the frequencies of the mixed modes. 

Regardless of which family of solutions is a better repre- 
sentation of KIC 11026764, the parameter values in Table [4] 



already yield precise determinations of some of the most inter- 
esting stellar properties — in particular, the asteroseismic age 
and radius. If we consider only the four models that were 
produced by the same modeling team with an identical fitting 
method (the ASTEC models: AA, AA', AB, AB'), we can cal- 
culate the mean value and an internal (statistical) uncertainty 
in isolation from external (systematic) errors arising from dif- 
ferences between the various codes and methods. The astero- 
seismic ages of the two families of solutions range from 5.87 
to 5.99 Gyr, with a mean value of t = 5.94 ±0.05 Gyr. The 
stellar radii of the two families range from 2.03 to 2.08 Rq, 
with a mean value of R = 2.05 ±0.03 Rq. The luminos- 
ity ranges from 3.45 to 3.80 Lq, with a mean value of L = 
3.56 ± 0. 14 Lq. These are unprecedented levels of precision 
for an isolated star, despite the fact that the stellar mass is still 
ambiguous at the 10% level. Of course, precision does not 
necessarily translate into accuracy, but we can evaluate the 
possible systematic errors on these determinations by looking 
at the distribution of parameter values for the entire sample of 
modeling results, not just those from ASTEC. 

Considering the full range of the best models (\ 2 < 10) in 
Table HI there is broad agreement on the value of the stellar 
radius. The low value of 2.03 Rq is from one of the ASTEC 
models considered above, while the high value of 2.09 Rq 
is from the Catania-GARSTEC code, leading to a systematic 
offset of ±[}q2 Rq compared to the mean ASTEC value. There 
is a slightly higher dispersion in the values of the asteroseis- 
mic age. Again considering only the best models, we find 
a full range for the age as low as 4.99 Gyr from Catania- 
GARSTEC up to 5.99 Gyr from ASTEC, for a systematic 
offset of ^ f 5 Gyr relative to the ASTEC models. The best 
models exhibit the highest dispersion in the values of the lu- 
minosity, with a low estimate of 3.45 Lq from ASTEC and a 
high value of 4.44 Lq from Catania-GARSTEC. These mod- 
els establish a systematic offset of ^'jf Lq compared to the 
ASTEC results. 

To evaluate the relative contribution of each observational 
constraint to the final parameter determinations, we can 
study the significance S-, using singular value decomposition 
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(SVD; see iBrown et all [19941) . An observable with a low 
value of Si has little influence on the solution, while a high 
value of St indicates an observable with greater impact. The 
significance of each observable in determining the parameters 
of KIC 1 1026764 for Model AA is shown in Figure H From 
left to right we show the spectroscopic constraints followed 
by the I = 0, 1 and 2 frequencies labeled with the reference 
radial order from Table Q] and each group of constraints is 
shown in a different shade. It is immediately clear that the 
I = 1 and I = 2 frequencies have more weight in determining 
the parameters. If we examine how the significance of each 
spectroscopic constraint changes when we combine them 
with modes of a given spherical degree, we can quantify the 
impact of each set of frequencies because the information 
content of the spectroscopic data does not change. The spec- 
troscopic constraints contribute more than 25% of the total 
significance when combined with the / = modes (13% from 
the effective temperature alone), indicating that these two sets 
of constraints contain redundant information. By contrast, 
the total significance of the spectroscopic constraints drops 
to 4% for the / = 2 modes, and 7% when combined with the 
/ = 1 modes — confirming that these frequencies contain more 
independent information than the I = modes, as expected 
from the short evolutionary timescale for mixed modes. 
Although the significance of the atmospheric parameters on 
the final solution appears to be small, we emphasize that 
accurate spectroscopic constraints are essential for narrowing 
down the initial parameter space. 

6. SUMMARY & DISCUSSION 

We have determined a precise asteroseismic age and radius 
for KIC 11026764. Although no planets have yet been de- 
tected around this star, similar techniques can be applied to 
exoplanet host stars to convert the relative planetary radius 
determined from transit photometry into an accurate abso- 
lute radius — and the precise age measurements for field stars 
can provide important constraints on the evolution of exo- 
planetary systems. By matching stellar models to the in- 
dividual oscillation frequencies, and in particular the I = 1 
mixed mode pattern, we determined an asteroseismic age and 
radius of t = 5.94 ± 0.05(stat)!° ;^(sys) Gyr and R = 2.05 ± 
0.03(stat)±QQ2( s y s ) Rq. This represents an order of magni- 
tude improvement in the age precision over pipeline results — 
which fit only the mean frequency separations — while achiev- 
ing comparable or slightly better precision on the radius (cf. 
Models J' and K' in TableHJl. The systematic uncertainties on 
the radius are almost negligible, while the model-dependence 
of the asteroseismic age yields impressive accuracy compare d 
to other age indicators for field stars (see ISoderblom 20101) . 
Whatever the limitations on absolute asteroseismic ages, stud- 
ies utilizing a single stellar evolution code can precisely de- 
termine the chronology of stellar and planetary systems. 

The I = 1 mixed modes in KIC 1 1026764, shifted from reg- 
ularity by avoided crossings, play a central role in constrain- 
ing the models. Bedding et al. (in prep.) have pointed out 
the utility of considering the frequencies of the avoided cross- 
ings themselves, since they reflect the g-mo de component of 
the ei genfunction that is trapped in the core (lAizenman et alJ 
119771) . The avoided crossing frequencies are revealed by the 
distortions in the / = 1 modes, which are visible as rising 
branches in Figure Q] For Model AA, marked by the verti- 
cal line, we see that the frequencies of the first four avoided 
crossings are Gi « 1270 /iHz, G2 ~ 920 /iHz, G3 ~ 710/iHz 



and G4 ~ 600 /iHz. Each of these avoided crossings pro- 
duces a characteristic feature in the echelle diagram that can 
be matched to the observations. Indeed, the observed power 
spectrum of KIC 11026764 (greyscale in Figure 0) shows a 
clear feature that matches G2 and another, slightly less clear, 
that matches G3. The avoided crossing at G\ that is predicted 
by the models lies outside the region of detected modes, but it 
is possible that additional data expected from the Kepler Mis- 
sion will confirm its existence. We also note a peak in the 
observed power spectrum at 586 /iHz (greyscale in Figure [3) 
that lies close to an / = 1 mixed mode in the models, and also 
at 766 /iHz near an / = 2 mixed mode. Again, additional data 
are needed for confirmation. 

It is interesting to ask whether all of the models discussed 
in this paper have the same avoided crossing identification. 
For example, are there any models that fit the observed fre- 
quencies but for which the avoided crossing at 920 /iHz cor- 
responds to Gi instead of G2? This would imply a different 
local minimum in parameter space, and a different location 
in the p-g diagram introduced by Bedding et al. Although 
this may be the case for some of the models in Table [3] from 
the initial search, all of the models listed in Table H] have the 
avoided crossing identification described above. 

The observed structure of the I = 1 ridge suggests relatively 
strong coupling between the oscillation modes. At frequen- 
cies above the observed range, the best models suggest that 
the unperturbed / = 1 ridge would align vertically near 40 /iHz 
in the echelle diagram (see Figure |3J- It is evident from Fig- 
ureQ]that at a given age the frequencies of numerous p-modes 
are affected by the rising g-mode frequency, and this mani- 
fests itself in the echelle diagram with several modes deviat- 
ing from the location o f the unperturbed ridge on either side 
of the avoided crossing (Deheuvels & Michel 2009). Stronger 
coupling suggests a smaller evanescent zone between the g- 
mode cavity in the core and the p-mode cavity in the enve- 
lope. Although the evanescent zone will be larger for I = 2 
modes, the strength of the coupling for the I = 1 modes raises 
the possibility that weakly mixed I = 2 modes — like those near 
766 /iHz in Model AA — may be observable in longer time se- 
ries data from continued observations by Kepler. 

It is encouraging that with so many oscillation frequencies 
observed, the impact of an incorrect mode identification ap- 
pears to be minimal. For example, the models suggest that 
the lowest frequency I = mode in TableQ]may actually be on 
the I = 2 ridge — or it could even be an / = 1 mixed mode pro- 
duced by the G4 avoided crossing. However, the influence of 
the other observational constraints is sufficient to prevent any 
serious bias in the resulting models. Even so, adopting either 
of these alternate identifications for the lowest frequency / = 
mode would cut the \ 2 of Model AB nearly in half. 

Despite a 10% ambiguity in the stellar mass, we have 
determined a luminosity for KIC 11026764 of L = 3.56 ± 
0.14(stat)^J jj(sys) L Q . With the radius so well determined 
from asteroseismology, differences of 200-300 K in the effec- 
tive temperatures of the models are largely responsible for the 
uncertainties in the luminosity. These differences are gen- 
erally correlated with the composition — hotter models at a 
given mass tend to have a higher helium mass fraction and 
lower metallicity, while cooler models tend to be relatively 
metal-rich. The adopted spectroscopic constraints fall in the 
middle of the range of temperatures and metallicities for the 
two families of models, and leave little room for substantial 
improvement. Perhaps the best chance for resolving the mass 
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ambiguity, aside from additional asteroseismic constraints, is 
a direct measurement of the luminosity. Although Kepler was 
not optimized for astrometry, it will eventually provide high- 
quality parallaxes dMonet et al.l2010f) . The resulting luminos- 
ity error is expected to be dominated by uncertainties in the 
bolometric correction (^0.02 mag) and the amount of inter- 
stellar reddening (^0.01 mag), although saturation from this 
bright target may present additional difficulties. This should 
lead to a luminosity precision near 3%, which would be suffi- 
cient to distinguish between our two families of solutions for 
KIC 11026764. 
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APPENDIX 
A. THE INFLUENCE OF ROTATION 

Rotation velocity can have a large influence on the frequencies of a given stellar model, since angular momentum transport can 
change the internal structure and evolution of the star. Thus, rotation velocity and angular momentum transport processes must 
be taken into account for accurate modeling. However, if the rotation velocity is not too large, the effect on the frequencies can be 
comparable to or smaller than the observational accuracy, in which case non-rotating models can be trusted. We have examined 
two angular momentum transport processes covering the extreme cases to quantify the influence of rotation on the modeling of 
KIC 11026764. 

The study was performed using the CESAM code (Morel & Lebreton 2008). The first-order effects of rotation on the equilib- 
rium models were considered by subtracting the spherically averaged contribution of the centrifugal acceleration to the gravity of 
the model, g e ff = g — A c (r), where g corresponds to the local gravity, and A c (r) represents the radial component of the centrifugal 
acceleration. This spherically averaged component of the centrifugal acce leration does not change the order o f the hydrostatic 
equilibrium equations. Such models are referred to as pseudo-rotatin g (seelSoufi et al.l fl998: Sua rez et al.l l2007). Since we have 
only a weak constraint on the rotation velocity of KIC 1 1026764 (see 83.21 i. we used an initial rotation velocity that leads to solar 
rotation at the solar age. Standard physical inputs were used, including the EFF equation of state. The opacity tables were taken 
from the OPAL package ([Rogers & Igle sias 1995), complemented at low temperatures (T < 10 4 K) with the tables provided by 
I Alexander & Ferguson! (fl994l) . The outer boundary conditions were determined by assuming a plane-parallel Eddington grey at- 
mosp here. The model metallicity (Z/X) is derived from the observed [Fe/H] value assuming (Z/X) Q = 0.0245 dGrevesse & Noels! 
1993), T pl = 0.235 and Z pl = for the primordial helium and heavy-element abundances, and AY / AZ = 2 for the enrichment ratio. 
The thermonuclear reactions incorporated the PP and CNO cycles with the NACRE coefficients. No microscopic diffusion was 
included in the calculation. 

We studied the following angular momentum transport processes: ( 1 ) Global Conservation of the angular momentum [GC, solid 
rigid rotation] : 51(f , r) = fl(t), and (2) Local Conservation of the angular momentum [LC, differential rotation] : d r 2 il/dt = 0. The 
GC and LC of angular momentum represent the two extreme cases in nature. The actual rotation profile of the star must fall 
between these two solutions. 

The theoretical frequencies were calculated using GraCo (|Moya & Garridol 120081) . We found that for the rotation velocity 
studied here, rigid rotation yields frequencies closer to the non-rotating case, with differences in the range [-0.01, -0.025] /jHz. 
Differential rotation yields larger differences in the range [-0.051, -0.045] /iHz. Considering the large frequency separations, 
rigid rotation leads to differences relative to the non-rotating case in the range [+0.0025, -0.0045] /iHz while differential rotation 
yields differences in the range [0, -0.004] pMz. All of these differences are much smaller than the observational errors, so 
rotation can safely be neglected in the frequency analysis of KIC 1 1026764. However, note that rotational mixing — which has 
not been considere d in this study — m ay lead to changes in the global and internal properties of the models even for slowly 
rotating stars (Eggenbergeret al. 2010l). 



B. THE IMPORTANCE OF ASTEROSEISMIC CONSTRAINTS 

We compared results from SEEK with and without the asteroseismic inputs to see how the large and small separations can help 
reduce the uncertainty on the inferred stellar properti es. We see that the uncertainty in the radius of 3% provided by SEEK is 
in line with what was expected from simulated data (Stello et al. 2009a). With R = 2. 10 ±0.06 Rq, the precision is a factor of 
five better than what we could get using only the available spectroscopic input (T e s, \ogg and [Fe/H]), which resulted in a 14.6% 
uncertainty (R = 2.06 ± 0.30 Rq). This dramatic improvement was obtained using the large and small separations as the only 
seismic constraints, as opposed to fitting individual mode frequencies. 

In Figure[5]we show 2D projections of the probability distributions provided by SEEK, which illustrate the correlations between 
the spectroscopic input parameters and the value of the inferred radius. To produce these figures we fixed all but one of the 
spectroscopic input parameters to the observed values. Hence, in Figure |5Ji the dark red regions are where we would most 
likely find the star if we did not know its temperature. The cross marks the observed value of the temperature with its measured 
uncertainty in the vertical direction while the horizontal direction marks the inferred value of R and its precision. We can see 
that changing the temperature within the error bars does not greatly affect the value of the radius. Reducing the uncertainty of 
the temperature would not affect the uncertainty of the radius, since the width of the correlation function is large compared to its 
slope. In Figure|5]5 we see that the value of logg is well determined without any spectroscopic measurement. The spectroscopic 
value of \ogg has no significant influence in constraining the radius. The constraints from the seismic measurements are simply 
much stronger than what is found from spectroscopy in this case. We found similar results for the determination of the mass. 
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Figure 5. Normalized probability distributions, showing the correlation of the radius with several spectroscopic constraints when they are excluded from the fit, 
including (a) the effective temperature T s {{, (b) the surface gravity logg and (c) the metallicity [Fe/H] relative to solar, where the gaps between the islands are 
caused by the discrete nature of the grid. Dark red indicates the maximum probability, which decreases linearly to white for zero probability. In each panel, the 
cross shows the observed values from the VWA method (see Table|2j and the final fit for R = 2.10 ± 0.06. We see that log j; is determined precisely before any 
spectroscopic knowledge of its value is included in the fit. 

Quite apparent in Figure[5j; are several local maxima due to the discrete nature of the grid. However, the resolution of the grid is 
sufficient to see the underlying correlation between the metallicity and the radius. From this figure we see that if the star had a 
metallicity of [Fe/H] ~ -0.6 it would have a radius of R ~ 1.91 7?©. We remark that the uncertainty in radius is mainly caused 
by the unknown value of the initial helium c ontent Y\ which has not been constrained by any o f the observables (for more details 
on the effect of Y| on stellar parameters, see Ouirion, Christensen-Dalsgaard & Arentoft 2010). 

Finally, in Figure |6]both the large and the small separations are left as free parameters. We note that the curved shape of the 
probability function is typical when the correlation between the large separation and the radius is plotted. In addition to showing 
that the large separation is the main constraint on the radius, Figure [6] illustrates that a star with a large radius, and hence smaller 
value of the large separation, will have a larger uncertainty on the inferred radius. The larger uncertainty is caused by the changing 
slope of the correlation, which tends to be flatter for larger radii. Figure|6]also shows that without the asteroseismic constraints, the 
uncertainty in the radius increases dramatically since the most probable region for the radius, in red, spans 1 .76 < R/Rq < 2.36. 

For the mass, SEEK demonstrates that the precision is only slightly improved by including the large an small separations in 
the fitting process. We find M = 1 .25 ±0.13 M© when only the spectroscopic input parameters are used and M = 1 .27 ± 0.09 M© 
when asteroseismic inputs are added. For the age, the results of SEEK without the asteroseismic inputs yield a relatively weak 
constraint of 40% with t = 5.03 ± 2.02 Gyr. This uncertainty is reduced by nearl y a factor of two when the large and small 
separa tions are used, which gives t = 4.26 ± 1 .22 Gyr. If we compare this result with Ou irion. Christensen-D alsgaard & Arentoft! 
(2010), 30% error on the age is unusually large for a star having well measured values of Av and 8v. It seems that 5v does not 
provide an additional constraint on the age of the star, since we were able to determine a very similar age t = 4.21 ± 1.25 Gyr 
when Sv was excluded from the fit. Bedding et al. (in prep.) find this to be generally true for subgiant stars. 

Finally, we stress that the mixing-length parameter a and the initial helium content F; are not constrained by our fit. For Yi, the 
output value is simply the central value used in SEEK's grid plus a la uncertainty extending toward the edges of the grid. For a, 
a value smaller than 1.8 is preferred, but not decisively. Our overall conclusion from SEEK is that the asymptotic asteroseismic 
inputs can be used to reduce significantly the uncertainty on the radius of KIC 11026764. The age and the mass precision can 
also be improved, though less substantially. 
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Figure 6. Similar to Figure|3] but for Ai^o when the small separation is also excluded from the fit. The point with error bars indicates the result when both the 
large and small separations are included as constraints. 



